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ABSTRACT 


The distribution of a measure of distance between a probability 
density function and its estimate is examined through Monte Carlo 
methods. The estimate of the density function is that proposed by 
Rosenblatt using sums of weight functions centered at the observed 
values of the random variables. The weight function used in all 
cases was triangular, but both uniform and Cauchy densities were 
tried for different sample sizes and bandwidths. The simulated dis- 
tributions appear in all cases to be close to Gamma distributions, 
but it has not been possible to relate the parameters to the popu- 


lation characteristics or to the window shape and bandwidth. 
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1, INTRODUCTION 


There are several recently proposed classes of empirical proba- 
bility density function [Refs. 1-5] all generally considered to be 
superior to the classical histogram estimates. The class considered 
in this paper is based on independent observations, 1.е., X), х„,..., 
x, are independent and identically distributed random variables with 
continuous unknown density function f(x). The method used to esti- 
mate f(x) is that proposed by Rosenblatt [4] and is as follows. 


Denote the estimator of f(x) by f (x). Then 





1 n x = X. 

= wJ —— 

3) 2 | b(n) | í 

where W(u) is a bounded integrable weight function with 


со 


/ W(u) du = 1 


-00 
and b(n) is a bandwidth that tends to zero as n > œ , but is such 


that o(b(n)) = Au i.e. it converges at a rate slower than 1/п. 


Thus we might have b(n) ~ 17/2. Note that all estimates of this 
form for given observations are themselves density functions; that 
is 


Е С) > 0, 


со 


/ f(x) dx = 1. 


«ОО 





Since the Х "5 are random variables, f(x) is a continuous parameter 
stochastic process. 

The local properties of such estimates satisfying relatively 2 
conditions can be shown to be biased [3]. Our object in this thesis 
is to investigate a global measure of how good Е (x) is as an esti- 
mate of f(x). This measure was proposed by Bickel and Rosenblatt [2]. 
The particular global measure ಹ್‌ for investigation is the 


‘Statistic: 


[£ (x) = £(x)1° | / 
В (а) -/ — а, 


whose value will vary — each realization of х,, 552; x Раме ТЕ 
is a statistic or function of the n random variables. 

Since an exact distribution of this statistic is not mathematically 
tractable, we have examined some representative cases by simulation. 
The conjecture made by Bickel and Rosenblatt [2] is that the 
statistic В(п) 15 asymptotically normally distributed regardless of 
the underlying distribution f(x). With this conjecture in mind, we 
hoped to establish a criterion that can be used to perform goodness 

of fit testing. 

A problem with this type of criterion is that, unlike the Kolmo- 
gorov-Smirnov or Cramer-von-Mises statistics for goodness of fit of 
an empirical distribution function to the unknown distribution, ß(n) 
is not distribution free and its distribution for finite n is unknown 
for any case. However, it is believed that ß(n) is fairly robust 
with rapid convergence to the asymptotic distribution; another object 


of the thesis is to examine this hypothesis. 





Going beyond asymptotic normality of f(n), the following result 
for a quadratic functional is also of some interest [2]. The func- 
tion a(x) is assumed to be a bounded piece-wise smooth integrable 


/ 


function. If b(n) = O(n? P) ав п”»%, 


Bun) > [nb(n) f (8, 69 - Є(х)17 а(хуйх - | ө сх)а(х)ах ЇГ н(а) 42] 


is asymptotically normally distributed with mean zero and variance 


21 19) (о) fo’ f(x)” dx 


(4) 


as n >°% , where W ls the 4-th convolution of W with itself. 





ІІ. SIMULATION 


The primary object of the simulation is to investigate the dis- 


tribution of the statistic В(п): 


[£ (x) - #0017 
В (п) -/ ತಾ = ах, 


Over a suitable range of integration. We have broken the simulation 


into the following cases: 


А. UNIFORM (0,1) CASE 
Since (287 0 <х < 1 
Е (х) = 


о, otherwise 


B(n) becomes 


1-Ь (п) 
В (п) -/ [£ ೧0-35 ах. 
Ь (п) 
We integrate from b(n) to 1-b(n) instead of from O to l, simply 
because we want to avoid the marked bias near O and 1. It can easily 


be shown that t, (x) is biased in the end regions, but not otherwise. 





Thus the mean and variance of Е (x) in this case 15: 


1 | 
а. х-у 
Elf (x)] = 575 Jute) dy | 
6 | 
х+Ь (п) 
1 |х-у | 
b (n) / Са b(n)? dy 
х-Ь (п) 
1 | х+Ь (п) — 
"ನಾರ "а ЧЭ рае] bin) Y 
x-b (n) х-Ь (1) х 
1 E ea 1 Фа)? 
вето З коро оо І: 


This is only true for b(n) < x < l-b(n). Also 








i п ү, X-X j) 
var Ыы 2, = | 


3-1 


цал 


~¥ (5) 
ಸಜ ಗರು) 
nb (n) 
1 


ato Ue ванн (Ля) „): 


೦ 


п ne 222 (п) 


The integration was performed in the simulation using Simpson's 
rule, dividing the interval of integration into 100 equally spaced 
sub-intervals. The results were found to be satisfactory in the 
sense that the value of the integral for 100 sub-intervals differs 


little from those with more sub-intervals. 





To evaluate f (x), use the following formula: 





х-Х,. 
£ ೧0 ಇಲ್ಲಿದೆ ul ОМ 
п пЬ (п) за Ь (п) 


where the weight function W(u) is triangular, і.е.: 


14, if [ч | < 1 
W(u) = | 
0 : otherwise 


Kr j=1, ..., n are independent and identically distributed random 
variables. In this first case, they aré uniform (0,1) random variables. 
Also b(n) is the bandwidth. We employed three bandwidths in the 
uniform (0,1) case, namely: 3/va, 1/уп, 1/n, One at a time. For each 
bandwidth, we used five different sample sizes. 1.е. 100, 200, 500, 
1000, 1500, and for each sample size, we used five different random 
number seeds, i.e. simulated five times. The reason for this is to 
have the ability to assess the variability of the simulation results 
from the sections. The random number seeds are taken from the imme- 
diate previous "1500 sample size" output, except the first random 
number seed which is arbitrary. So, there are five sets of random 


numbers which are entirely different from one another. The following 


table should help to clarify the scheme: 


10 ' 


> В Р 








TABLE 1 


Bandwidth Sample size, n No. replications # Sections 





2000 
2000 
2000 
2000 
2000 


ಟಾಟ (л (л (л 


2000 
2000 
2000 
2000 
2000 


ಟಾಟ (л (л (л 





100 2000 
1 200 2000 
п 500 2000 
1000 2000 
1500 2000 


сл (л (л (л (л 


For each sample size, n, repeat five times using different random 
number eo. 

For the uniform (0,1) case, there are 75 outputs altogether. In 
- each single simulation (one computer output), we evaluated 400 ß(n)'s, 
say В. (п), j=l, «ee, 400. The output includes the 211221 

l. A bar histogram of the 400 В(п) "5. 

2, SNAP/IEDA graphs: 

a. A refined solid histogram. 


b. Graph of empirical log-survivor function, &£n(1-j/(N+1)) vs 


В (у) (п), where В (n) is an ordered value. 


3) 


ne 


W“ 


241 








с. Empirical CDF graph i/(ntl) vs Bij) (СО, where jy (n) 
denotes the ordered observed В. (п) "8. 

4. Normal probability distribution with estimated mean and 
variance. Norm vs e.c.d.f. of В, (п), where Norm is computed as 
follows: Let 


У В. (п) 
A bi 


n 


and 82 (В(п)) = sample variance of the 8. (а). Then 


| | ЖЫ - р 
Мои ал: | S(B(n)) | | 


e. Normal probability plot Norm vs i/(N+1). 
A sample output is included in the "Computer Output" section 
(page 30 - 35). 

For each bandwidth and each sample size, there are five computer 
outputs, or 5 x 400 = 2000 Р. (п) "8. With this 2000 sample size, we 
constructed a histogram again. Also we used the Rosenblatt smooth- 
ing eechniaue to get a smooth ಗಜಾ rima estimate for fín). 
After a number of trials, we found that the bandwidth needed in the 
Rosenblatt estimator has very much to do with the smoothing and how 
"good" the density looks. Fortunately it is robust to a certain 
extent. This graph can give us a better picture of the true distri- 
bution than the classical histogram. The last graph of this complete 


output for sample size 2000 is the empirical log survival function; 


£n(1 - i/(N+1)) vs В (п). 


(3) 
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There are 15 "2000 sample size" computer outputs in the uniform (0,1) 

сазе. (Since each bandwidth has five "2000 sample size" outputs.) 
The complete computer outputs for each bandwidth and each sample 

size in the uniform (0,1) case are included in the "Computer Output" 


section (page 42 - 86). / 


B. THE CAUCHY DISTRIBUTION 
The next case investigated is the Cauchy distribution, again with 


a triangular weight function and with the range of integration (-3,3): 


+3 2 
Ге, (x) - £(x)] 
B(n) = Fm ах : 
-3 
The interval of integration comprises 80% of the probability mass for 
this distribution. The integration is carried out, in this case, 
using Simpson's rule with 600 sub-intervals. This grid was used after 
examining the 100, 300, 600, 900 sub-interval cases to ensure that a 
fine enough grid is used in the numerical integration. 
The Cauchy density function is 


f (x) =, 


2 
T(1 + х ) 
The density function estimator is the same as in the uniform case: 


Х-Х, 


ಷ್‌ Е нів ° 


The weight function W(u) is again triangular. Thus the only differ- 


f, (x) = 





ence in this formula from the previous case is that Kae j=l; 25.2; N 
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are i.,i,d.Cauchy random variables. Note also that for finite n the 
estimator will have a bias component, this bias component usually 
decreasing, for given n,as bandwidth decreases. In that case though 
the (pointwise) variance of f (x) increases and possibly also the 
variance of ß(n). Note that ß(n) and its distribution is a function 
of both sample size and bandwidth. 

The bandwidth b(n)'s used here are 1/уһ, 3/10; 20/Vn, the last 
one representing a case in which bias in the estimate f (х) plays a 
major role in a distribution of ß(n). 

The rest of the simulation scheme is exactly the same as for the 
uniform (0,1) case. 

A sample output for the 400 В; (п) case is included in the "Com- 
puter Output" section (page 36 - 41). Also the complete computer 
outputs (2000 replications) for each bandwidth and each sample size 
in the Cauchy case are included in the "Computer Output" section 


(page 87 -131). 


C. HAND PLOT GRAPHS 

Besides the computer work, there are six hand plot graphs. These 
are the ß(n) quantiles versus the standard normal quantiles, They 
are included in the Appendix, to illustrate the rate of convergence 


of the distribution of 8(n) to its asymptotic normal distribution. 


14 | \ 





111. RESULTS 


The asymptotic result, conjectured by Rosenblatt, is that by 
using a triangular weight function when estimating the uniform 


density function on (0,1) 


1-Ь (п) 
-1/2 2 2 
{b(n) } nb (n) |£ (x)-1|%4x-(1-2b(n)) | Wu) “du 


Ь (п) 


і 


is asymptotically normally distributed with mean zero and variance 


ow ©) (0) (1-2b (n) ) 


2/9 


as n > œ if nb(n) > © and b(n) = on ). Here n is the sample 


size and b(n) the bandwidth. For the triangular weight function 


со 


wie) du = = 


and и (о) the 4th convolution of W with itself at zero, is 302/630. 
From the above expressions, we get: 
1-Б (п) і 
2 Е 2 (1-2b(n)) 
E(B(n)) = ГЈ |“ (х)-1| “|. иж 
Ь (п) 
1-Ь (п) (4) 
eve 1Е (х)-112ах| . 28 (0) (i-2b(n)) 
n 2 
n b(n) 
Ь (п) 


A comparison of the simulated values with that of the conjectured 
is tabulated below for the different cases investigated when f(x) is 


a uniform density: 





TABLE II 


Comparison of estimated mean values and 
asymptotic mean values of В(п) for dif- 
і ferent bandwidths and sample sizes. 


Е (В (п) ) / (1-2b (n) ) 
Б(п)=3/Уп Е(В(п)) Conjectured Computer output 


.3000 »0089 .0222 . 0127 
.2121 «0090 „0157 • 0109 
„1342 „0073 „0099 • 0075 
• 0949 • 0057 „0070 „0058 
„0775 „0048 „0057 • 0051 


Е (В (п))/ (1-2ь (п)) 
Ь (п) =1/Уп Е (8 (п) ) Conjectured Computer output 





TABLE III 


Comparison of estimated standard deviation 
values and asymptotic standard deviation 
values of ß(n) for different bandwidths and 
sample sizes. 


| с (В (п)) / (1-25(п)) 
Ь(п)-3/Уп С (В (п) ) Conjectured | Computer output 


• 3000 • 0113 
„2121 • 0081 
• 1342 • 0046 
. 0949 „0029 
„0775 „0022 


0(8(n))/(1-2b(n) ) 
Ы(п)-1/Уп С (В (п) ) Conjectured | Computer output 
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а Тог small bandwidth the agreement between the asymptotic 
and simulated variances is very good even for small n(n=100). The 
same is true for the expected value, although convergence is slower 
than for 0 and again slower for large bandwidth. 

The simulation reveals that the distribution of the 8(n) statistic 
converges to normality very slowly (cf. the six ೬ in the 
Appendix). For sample sizes 100, 200, 500, 1000, 1500, the histograms 
are all skewed to the right; the empirical log-survivor function plots 
are concave downward; and none of the normal probability plots is 
approximately straight lines (refer to "Computer Output" section). 

It is very suggestive that the 8(n) statistic is approximately Gamma 


(0,k) distributed where the density of the Gamma distribution is: 


к. = 
2. 0812-1222 
(К-1)! 
and the mean and variance are: 
. 2 
E(X) - 0; Var(X) - ко“ . 


Therefore, for each underlying distribution, each bandwidth, and 
each sample size, we used the 2000 evaluated f(n)'s to estimate the 


0's and k's, and tabulate them as follows: 


17 \ 








TABLE IV 


Estimated parameters for fitted Gamma distribution for 
B(n) in the case where the densities are uniform (0,1) 
and Cauchy, with a triangular weight function. 


(1) Bandwidth = 1/n Underlying distribution: Uniform (0,1) 


Sample sizes and estimated 8's and k's 


| 100 200 1500 (std. dev.) 


.6540 | .6638| .6612| .6621 „6645 ж  .0045 





к | 39.48| 38.46| 32.95| 30.86 31.06 + 2.91 


(2) Bandwidth = 1/Уп Underlying distribution: Uniform (0,1) 





Sample sizes and estimated 8's and k's 


| р 100 | 200 | 500 1000 1500 (std. dev.) 
8 | .0565] .0420| .0279| .0201|  .o166 + .0003 | 


3.62 5.56 8.63! 12.61 16.92 + 1.45 


Seren sizes and — 8's and k's 









Underlying distribution: Cauchy 









1500 (std. dev.) 
» ој. СК ee 
К 21.02 31.49 59.34 79.19 100,83 + 4.37 

(3) Bandwidth = 3/vn Underlying distribution: Uniform (0,1) 


ple sizes sane estimates 0's and k's 
100 1000 1500 (std. dev.) 


A .0123| .0107| .0079| .0060 .0051 + .0001 
1.59 2.43 3.65 5.03 + 0.37 


Sample sizes and estimated 0'ѕ and k's 


Т 









Underlying distribution: Cauchy 


1500 (std. dev.) 


6 | ‚1247 • 0911 • 0576 • 0412 • 0340 • 0004 


К 8.62 12.30 20.32 29.14 33.98 1,55 





(4) Bandwidth = 20/Уп Underlying distribution: Cauchy 


le sizes and estimated = ane k's 
100 (std.dev.) 200 1000; 1500 (std. dev.) 


Н .0558 .0008 „0262 -0110|.00661.0052 • 0001 
6.71 + 0.78 4.101 5.24 + 0.16 


where 8 = mean , 





Sam 







1/(estimated coefficient of МАРТЫ ТО)" 


^ 
e 
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The standard deviation estimates for the parameter estimates are 
obtained from the five identical simulations using different random 
number seeds, each generating a sample size of 400, as detailed in 


Table I. 


19 





IV. CONCLUSIONS 


Although the conjecture proposed by Rosenblatt is that the dis- 
tribution of the statistic В(п) is asymptotically normal, this is 
certainly not the case for moderate sample sizes such as those used 
in the simulations. Based on this fact, the conjecture has little 
practical value, but the fact that all cases examined give rise to 
an approximate Gamma distribution gives some hope for stability in 
the results. 

However, the method proposed by Rosenblatt for estimating density 
functions is no doubt plausible as can be seen from the Rosenblatt 
type estimates used in this study. More effort in this area is 
needed, particularly to relate the Gamma distribution parameters to 
underlying factors. Since we realize that for moderate sample sizes, 
the statistic 8 (1) is not distribution free, and the bandwith Ь (п) 

(a function of sample size) plays a crucial role in estimating the 
density аа and also the shape of the weight function W(u) is 
a major factor for evaluating f (x), the future work may be to find 
the best combination of the above three factors by simulating the 
distribution of the statistic fín) for different underlying distri- 
butions, bandwidths, sample sizes and weight functions and tabulate 
the resulting quantiles. Then we may be able to put the statistic 


B(n) to actual use. 


20 | п. 





APPENDIX 


Included in this Appendix are six plots of sample size versus 
estimated quantiles of the statistic f(n), normalized by subtract- 
ing the estimated mean and dividing by the estimated standard 
deviation. Also plotted as straight lines are the corresponding 
quantiles of the unit normal distribution. 

The 17 ог extracted from each 2000-8(n) sample, are 
0.010, 0.020, 0.025, 0.05, 0.10, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70, 
0.80, 0.90, 0.95, 0.975, 0.980, 0.990. 


Other data for each plot are self-explanatory. 





= 
Ea — T E 


АЕ 
— — СО — ТЕ 
ЕЕЕ ваја “Врэг Т | 
—— 


En 
— ШИН 
РРР ЕНЕ ЕНЕ = 


FERE ЕЕ ады, 
= с анара эч HE TER 


фа 15 — 









ЧГ: 

і 

і 1 | 1.9 9 

ее Бе —— 
--0. 980 








~ 


— — 


шин Ч 
Е 


— rr — Do 80: 


ү Ж 
EEE —— = 
— — — 0.80 

— — — NN TT | 

| He Т Н 


— 






a] 


=> je T] 
= —— O le ШИ УВО | et ಜುಂ Ш ај 
cee. ent | 
e< AAA ka. „|“ ಸಮಾಜ. | эн | MIN | 
-1 | 4 —— — i М ; . 
АНЕ ше Шетен CE ES | чии | 
| | š — — Ша 
| : 
ವದು | ЕСЕ Еј 11 ы] 
೬1. ಸ್‌ —i po | SE 
— —— Осо: с соон [жыгыша сє ET 


— — ај 





о 


Estimated (n) quantiles 5 Н(О,1) quantiles 





шо» > © Sample size, n 


1000 


да, 





Estimated (n) quantiles & N(0,1) quantiles 


4 
= 


КЕГЕ ааа А У АНЕ 





| . 
a 
ji 
| 





É ' 
беа 
— 


НЕЕ 
— ЕТ 








! | 
ЕЕ 
Без! ರ್‌! 


! | “+. 4 2 1 : ia ва | — эн Ь pw ы ще: mo гэр A ж - = —— — т р-а T ದ 
в — m HIT 


—— — = 


- | я 
ы 
, 
а 
! 






шиг Ju 


54 — — г i 
— — — — 


a і 
а 
-i m =- — com ess s... ಇರಾ — — 
а . 
| 5 I 
y t 





к А. — — — —— | 


И 


. 





-2 [s= ерес ааа 
2 310 


-3 


4 
. 
| i 
q | š | ! 
„а ми as . = ed жа = -ha “hy ——— 


, 
' 
, 
t 


Е 1111 16 қаң 


= —» 
3 
, 





20 
24 


о 
с А Sample size, n 


43 





ae 


ЕЕ 
— 
JN 


Hr 
НЕ ЕН 
КЕ eee НН 
ы КЫКЕ "ЕЕЕ 


шэн їг її 
eso dea ЕНЕЙ [ಇ 
| А ОДА Р ҮГ 
ಕಜ — ಶಲ 
AE EEE жы зы == s ш _0,99 


е ене ЕН НО ПЕН ЕН А 
=> , ы | i 0.98 


— — — a 7 — I  . 0.975 


— — — ВЕ 


ಪಾಪು ЕЗДІҢ 50 


| мета а [oles БГ је | а 
| EEE Н; 


| азы ш у ою ою 2 7 | 


1 
гүү лүү — 
L i; je i `: ESS 


ка ы | | 
| = pepa ay оа 
ЗЭЭ с” — 


НЕ — А 


=== re БЕ ~! 
| | 5 
2-1 — ----- қаша ° +» - — т DEY ЩЕ. 
и т | 
| — на ba 
! = j 
1! — = sc ~ ee jet 
| E — Er 


сее Ер i 
Ы ТІРІНІҢ | 


О ೦ 
в © О © Sample size, n 

















Estimated in) quantiles & N(0,1) quantiles 












ВЕБ 





== 
о 
О 
О 
ен 


1500 


24 | | 





Estimated B(n) quantiles € N(0,1) quantiles 


СЗ НЕ ВЕ НИ ЕТ Тт И ВС СБ 
Вен Ши 11 11111111 11:01 11: e 
ее Нее 
| | | | ||] | чраћејио! ада 


ШОШ ШҮ he 
EEE aa EHH 





і ка 
ШЕНЧЕНЕНЕЕНЕНЕНЕЕН 
—— * eae e a 


еі разсее 515 


= 





— 0.980 


4 
! 


цэв тын сак йг Р? ГЭЭ үнэ a == 5 95 | 


— — 


Пэн Ж 


E E на E 
— аты. ша а 
0 


ж а —FR — — -8 
Peer | — — — 
a EG ера је Геја] —p. | 





Ko EE Et; — — — г 4 60. 
Ёст ವರಾ ಹ EEE на 

ЕЕ ЕН — — AAA CE “lie 20 | 
Eee ев 
ЕЕ г жо ат оо В aa a моја 
— ರಸ а 
EET к=н акш ш УТД 


Ең — — ог түг 


| 


Об 
E ЕЕС СТАНЕТЕ "а 


eS эээ ———=sk= ES 
5.029 


— ЛЕ EEE AEE: 
== SE —— 

m 4 
Бү — — 
ЈЕ Н ЛҮҮ 


о | 
с © S O Sample size, n 
- ಇ 


Гүээ CTE Ж 
ЕТЕШ ЕЕ НЕЕ 75 
TA 











5 
1000 
1500 


AS 








Р Гр Ер ЫЕ 
МИН 1 1 1 Пен 
| PEP ಬಯ ಜಬ ಚಮಚ! 
БТ 
ಆಗತಾ ಸರ ಸೂರ್ಯಾ ತು ಸಮ 
Set ІІМ" ET EEE EERE TE 
Eee EEE ат 
АЕ ee E EERE EE 
ES ಪೂರಾ ETE 
Багт A ae 
вена Fea TE AEE TEE 
= ECE FEE HEHEHE 
AE TTT TEE 
LTTE г 
ME 12001201 008 
el [kl јин ај ——— — FE 
Кн” ЕЗ а ಸ. 
as Ко 
— 


ЕЕ ЕР el 















Estimated (n) quantiles & N(0,1) quantiles 








2] eS и 


шшш ср ee Se ee ee шы ee | Ээ ee ee ee 


SEE H a 


p — a. R 


| ಸ ರರು z 
шиг EE TEE 


— de dest=} | ЈЕ | ај ЕТ ЈЕ ЈЕ ЈЕ | — — — — 
== = === 





== ще шиг. ло 
23 Сү тэ түү 51) 
- + 





A ганан ыкташы "алышына ——  — = 


ГГц ಸಃ 0.20! 
ШЕ; ані. 





— — — — I 

Ec к L — 

panel 
Eur “Piti a ry: 








8, 6 





+ wists ease . 
на | і и 
55 at t. 
AE НЬ ЧЕ 
+ NS . 





+ 





ыы ЕЕ Ili: 
нг 6 A 1 = 
= үе нин АЈ ЕНБ ИН | 

ಇಹ ೫3 Е АЕ те е аа Та ы ГЕ 
EE А А eee 
208 кке ЕН а ре š 


pen La се| токо сев 














Ss E EFI TE РЕ Е ЕЕЕ 





mtb - б ОЖ ЕЕ | ао 
ЯЕНЕНЕНЕНЕНЕР шэг ээ AO 
22. анамнан таны ыы 








| = БЕН E 
Bee ee 
ERE TTT ETE ET 
ЕТ шерт = з L 
Be ee me РЕ Б Е ЕЕ | 
: ШОГ: сг ао ಬುಸು 0.98 
2 Se en ee | < a 


A ја Teel lee 
Y —— — — — — — Го 4 
ӨГ ರ್ಸ್‌ шээс оо НИ. 
А ай __н A A AA чи 
Вер РАЈЕ 1 АРА ЕТ 


1 AL 5 ಎ : : 
— — — — ಹರರ a е! ~) 


ЕЕ — не 
БЕ ЕЕ Ер аер, 
ee Er: К, л PS 


ರ 
р ыс: ol Ме at al и: A 
ಹ BEE >< == 4% 


























Е ан шагаа ОШ. аламы! у Sy o A | у 0.40 
== АЕР М рар алата | 


деды а-а ыыы ааа а Р 


Estimated B(n) quantiles & N(0,1) quantiles 
— 





Ширээгээ л ың. 

O ಎಷು ಮಾವ 

Ве [| ПА Р С a ЕА | 

a гэ Ээ ЭЭ Эв СЭ Р 02: 
ZA E. A IA | 


EEE EET 
EEE EE 





тя ве e. 


.. El 
Es — — 4 -4 --- Ве / | - - --- — 
> ие за . . .. “ ae ve . . Е 1 I 4 | Қ 
.... 4-.%.... eorthesarfro. Poo» .. .. а» ооо руа оао фа .. , , 
Е С ...зйЁй. "ре .. .. a+ «Р. for fo... ” е ° 






о о о 
= о о O 
еч N ಓಗಿ 


Sample size, n 


1000 
so 


9.7 





COMPUTER OUTPUT 


There are altogether 102 sheets of graphs. All of them are self- 


explanatory, except that the letter "X" in the graphs represents the 


statistic В(п). 


The following is a list of index of the graphs. 


36 - 


41 


Description 


sample computer output for the case: 
underlying distribution-uniform (0,1); 
underlying distribution sample size n=1500; 
bandwidth 3/vn; 

statistic B(n) sample size m=400. 


sample computer output for the case: 
underlying distribution-Cauchy; 
underlying distribution sample size n=100; 
bandwidth 1/Уп; 

statistic f(n) sample size m=400, 


The rest of the computer graphs in this section are all for statistic 


B(n) sample size m=2000. The page index is as follows... 


Page 
42 - 
45 - 


48 - 


54 - 
57 - 


60 - 


44 
47 
50 
53 
56 
59 


62 


Underlying 

Underlying distribution 

distribution sample size Bandwidth 
0(0,1) 100 З/Уп 
0(0,1) 200 З/Уп 
0(0,1) 500 З/Уп 
0(0,1) 1000 3/vn 
0(0,1) 1500 з/п 
0(0,1) | 100 1/Ул 
0(0,1) 200 3/1೧ 
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Underlying 


Underlying distribution 

Page distribution sample size Bandwidth 
63 - 65 0(0,1) 500 1/уп 
66 - 68 0(0,1) 1000 1/Уп 
69 - 71 0(0,1). 1500 1/Уа 
72 - 74 0(0,1) 100 1/п 
75 - 77 0(0,1) 200 1/n 
78 - 80 0(0,1) 500 1/n 
81 - 83 0(0,1) 1000 1/n 
84 - 86 Bro 1500 1/n 
87 - 89 Cauchy 100 1/Ул 
90 - 92. Сачспу 200 1/Уа 
93 - 95 Cauchy 500 1/Уп 
96 - 98 Сачсһу 1000 1/Уп 
99 - 101 Cauchy 1500 1/vn 
102- 104 Cauchy 100 3/1೧. 
105- 107 Cauchy 200 3/vn 
108- 110 Cauchy 500 3/5 
111- 113 j Cauchy 1000 3/1೧ 
114- 116 Cauchy 1500 3/1೧ 
117- 119 Cauchy 100 20/Уп 
120- 122 Cauchy 200 20/Уп 
123- 125 Cauchy 500 20/1೧ 
126- 128 Cauchy 1000 20/Уһ 
129- 131 Сачсһу 1500 20/Уп 
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